A comparison between several correlated 
stochastic volatility models 

Josep Perello^'*, Jaume Masoliver^, and Napoleon Anento^'^ 

^Departament de Fisica Fonamental, Universitat de Barcelona, Diagonal, 647, 

E-08028 Barcelona, Spain 

^Gaesco Bolsa, SVB, S.A., Diagonal, 4^9, 08036'Barcelona, Spain 
Abstract 

We compare the most common SV models such as the Ornstein-Uhlenbeck (OU), 
the Heston and the exponential OU (expOU) models. We try to decide which is 
the most appropriate one by studying their volatility autocorrelation and leverage 
effect, and thus outline the limitations of each model. We add empirical research on 
market indices confirming the universality of the leverage and volatility correlations. 
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The multiplicative diffusion process is considered to be the most popular model 
in finance. We thus have that the log-price change without average describes 
the diffusive process X{t) = cr^i(t) (in the Ito sense), where a is the volatility 
(assumed to be constant) and is Brownian noise with {^i(t)^i{t')) = 6(t—t'). 
However, this model is unable to capture most of the statistical properties of 
real markets, those called stylized facts. 

The stochastic volatility (SV) models are a natural way out to avoid the 
inconsistencies of the log-Brownian model by still assuming that X follows a 
diffusion process but now with a random a. However, the dynamics of a has 
not been definitively attached to any specific SV model. In this sense, we want 
to discern between the most common SV models. Among others we have: 

(a) The Ornstein-Uhlenbeck (OU) SV model [1,2] 

t 

&{t) = -a{a - m) + ^^(t), (T{t) = m + k J e-''^'-''^^2{t')dt' . (1) 
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(b) The Heston SV model [3,4] (with V = a^) 



V(t) = -a{V -m^) + kVvUt), (2) 

t 

V{t)=m^ + k J e-'^'^'-''^ ^V{t>)^2{t')dt' . (3) 

— oo 

(c) And, finally, the exponential Ornstein-Uhlenbeck (expOU) SV model [5] 
(with Y = ln(cr/m)) 

Y{t) = -aY + -G(t), Y{t) =m + - [ e-'^^'-^'^^ityt' . (4) 
m m J 

— oo 

Note that the Heston model does not provide a closed equation for the volatil- 
ity. This can be a serious drawback in case we want to deal with analytical ex- 
pressions. In addition, we suppose that noises are correlated (i.e., {^i{t)^2{t')) = 
p6{t — t') with — 1 < p < 1). This allows us to explain some stylized facts. 

The leverage effect [6] is measured by: L(r) = {dX{t + TfdXit))/ {dX{tf). 
For the models above we have 



Lou(r) = 2p 



z/V2a (1 + u-^e 
(1 + z/^)^m 



e-"^ LH(r)~4p^e-"^ (5) 



i^expou (r) = ^e^V2 (2i.2e— - ar) , (6) 

for r > and L[t^ = for negative r. It seems impossible to derive an 
analytical expression for the Heston leverage because V{t) does not have a 
closed expression (cf. Eq. (3)). The leverage given in Eq. (5) is exact only when 
= k"^ /2a but Ref. [4] shows that this is not true for the Dow Jones. We will 
have thus to simulate the process to compute the leverage effect in each specific 
Heston case. Eqs. (5)-(6) are plotted and compared with data in Fig. 1. The 
data is too noisy to assert which is the most appropriate model. In any case, 
the three models have the desired exponential decay with a characteristic time 
of the order of few days. We also note that, without the correlation coefficient 
p, there is no leverage and data confirms that p < [7]. 

The volatility autocorrelation is C(r) = [{dX{tfdX{t+Tf)-{dX{tf){dX{t+ 
r)^)]/Var[(iX(t)^]. In this case, we can exactly derive this correlation for each 
model. We have 



Cou(r) = 



4z/2(2 + i/2) + i ^ 



CH(r) = i^^e— , 



(7) 
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Fig. 1. The leverage effect for several daily price indices. We also add the leverage 
function L{t) fit for the different SV models. The S&P 500 is much more liquid. 

, _ exp {4u^e-n - 1 
^cxpOvyT) — — 7-r^^ — ■ 



3exp(4i^2) _ I 



The volatility autocorrelation observed in empirical data shows at least two 
different time scales (see Fig. 2 and Ref. [8]). The shortest one coincides with 
the leverage correlation time while the second scale is of the order of years, 
around 10 times bigger than the leverage time. The OU and Heston models 
are not capable of reproducing the slowest time decay but the expOU model 
has a nontrivial decay with the two time scales: 



CexpOu(T) 



(ar > 1), CcxpOu(t") ~ e 



-(fc/m)2 



(ar<l), (9) 



where the faster decay is given by the diffusive coefficient k/m and the slow 
decay of the volatility is provided by the reversion coefficient a (cf. Eq. (4)). 

Therefore, we conclude that the expOU model appears to be as more realistic. 
Further investigations on this model are required to confirm this assertion. 
However, an alternative choice is to sophisticate the OU and Heston models. 
A second time scale can be generated by including a third SDK for m. In a 
previous work [8], we have studied this possibility for the OU model with some 
success by taking: m(t) = —amim — 0) + km^^it). A similar analysis can be 
performed for the Heston model. 
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Fig. 2. The volatility correlation for several daily indices. We add two exponential 
fits one with a short range (similar to the leverage time scale) and a second one for 
the larger time lags. 
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